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CONTROLLER SYNTHESIS FOR ROBUST INVARIANCE OF POLYNOMIAL 
DYNAMICAL SYSTEMS USING LINEAR PROGRAMMING 

MOHAMED AMIN BEN SASSI AND ANTOINE GIRARD 



Abstract. In this paper, we consider a control synthesis problem for a class of polynomial dynam- 
ical systems subject to bounded disturbances and with input constraints. More precisely, we aim at 
synthesizing at the same time a controller and an invariant set for the controlled system under all ad- 
missible disturbances. We propose a computational method to solve this problem. Given a candidate 
polyhedral invariant, we show that controller synthesis can be formulated as an optimization problem 
involving polynomial cost functions over bounded polytopes for which effective linear programming 

r) ■ relaxations can be obtained. Then, we propose an iterative approach to compute the controller and 

the polyhedral invariant at once. Each iteration of the approach mainly consists in solving two linear 
programs (one for the controller and one for the invariant) and is thus computationally tractable. 

,S^ • Finally, we show with several examples the usefulness of our method in applications. 



1. Introduction 



The design of nonlinear systems remains a challenging problem in control science. In the past decade, 

building on spectacular breakthroughs in optimization over polynomial functions [LasOU IPar03| . 
QQ I several computational methods have been developed for synthesizing controllers for polynomial dy- 

VO ■ namical systems [PPR041 ILHPTOS] . These approaches have shown successful for several synthesis 

problems such as stabilization or optimal control in which Lyapunov functions and cost functions 
C^ '_ can be represented or approximated by polynomials. However, these approaches are not suitable 

for some other problems such as those involving polynomial dynamical systems with constraints on 

states and inputs, and subject to bounded disturbances. 

In this paper, we consider a control synthesis problem for this class of systems. More precisely, given a 
k><( I polynomial dynamical system with input constraints and bounded disturbances, given a set of initial 

;_j ■ states P and a set of safe states P, we aim at synthesizing a controller satisfying the input constraints 

and such that trajectories starting in P remain in P for all possible disturbances. This problem can 
be solved by computing jointly the controller and an invariant set for the controlled system which 
contains P and is included in P (see e.g. |Bla99j ). We propose a computational method to solve this 
problem. We use parameterized template expressions for the controller and the invariant. Given 
a candidate polyhedral invariant, we show that controller synthesis can be formulated as an opti- 
mization problem involving polynomial objective functions over bounded polytopes. Recently, using 
various tools such as the blossoming principle |Ram89j for polynomials, multi-affine functions |BH06j 
and Lagrangian duality, it has been shown how effective linear programming relaxations can be ob- 
tained for such optimization problems [BGlOj . We then propose an iterative approach to compute 
jointly a controller and a polyhedral invariant. Each iteration of the approach mainly consists in 
solving two linear programs and is thus computationally tractable. Finally, we show applications of 
our approach to several examples. 
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2. Problem Formulation 

In this work, we consider a nonlinear affine control system subject to input constraints and bounded 
disturbances: 

(2.1) ±{t) = f{x{t),d{t)) + g{x{t), d{t))u{t), d{t) € D , u{t) € U 

where x{t) £ Rx ^ K" denotes the state of the system, d{t) £ D CI M™ is an external disturbance and 
u{t) € f/ C RP is the control input. We assume that the vector field / : M" xR"^ — ;• R" and the control 
matrix g : R" xR™ — >■ R'"'^p\ defining the dynamics of the system, are multivariate polynomial maps. 
We also assume that the set of states is a bounded rectangular domain: Rx = [a;i,i!ci] x • • • x [xn, x^] 
with Xfc < Xfc for all k € {1, . . . , n}; and that the set of disturbances D and the set of inputs U are 
convex compact polytopes: 

D = {de R""! aD,k ■ d < Po^k, V/c G /Cd} and U = {u (^ W\ au,k ■ u < Pu,k, VA; e ICu} 

where a/j./fc G R™", Pd^^ £ R, au^k £ ^^, Pu,k £ ^, ^D and Ku are finite sets of indices. We will denote 
by Rd = [di_-,di\ x • • • x [dmidm\ the interval hull of polytope D, that is the smallest rectangular 
domain containing D; and by Vx = {a;i,xi} x • • • x {x„,x^} and Vd = {diidi} x • • • x {dm-,dm} 
the set of vertices of Rx and Rd- The present work deals with controller synthesis for a notion of 
invar iance defined as follows: 

Definition 2.1. Consider a set of states P C Rx and a controller h : Rx — > C/, the controlled 
system 

(2.2) ±{t) = f{x{t),d{t)) + g{x{t), d{t))h{x{t)), d{t) G D, 

is said to be P-invariant if all trajectories with x(0) G P satisfy x{t) G P for all i > 0. 

Let us remark that this is a notion of robust invariance since it has to hold for all possible distur- 
bances. Let P_Q P Q Rx be convex compact polytopes. In this paper, we consider the problem of 
synthesizing a controller h for system (|2.ip such that all controlled trajectories starting in P remain 
in P forever. This can be seen as a safety property where P is the set of initial states and P is 
the set of safe states. The problem can be solved synthesizing jointly a controller and a polyhedral 
invariant P C Rx containing P and included in P: 

Problem 2.2. Synthesize a controller h : Rx -^ U and a convex compact polytope P such that 
P_(l P (1 P and the controlled system (12. 2p is P-invariant. 



In the following, we describe an approach to solve this problem. To restrict the search space, 
we shall use parameterized template expressions for the controller h and the invariant P. Firstly, 
we will impose the orientation of the facets of polytope P by choosing normal vectors in the set 
{7^; G R"| k G fCx} where Kx is a finite set of indices. Then, polytope P can be written under the 
form 

P = {x G R"| 7fc • X < %, V/c G Kx} 
where the vector rj G RI'^-^I, to be determined, specifies the position of the facets. The facets of 
P are denoted by F^ for k G K-x-, where P^ = {x G R"! Jk ' ^ = Vk and 7i • x < rji, Mi G Kx \ {k}} ■ 
For simplicity, we will assume that the polytopes P and P are of the form: P = {x G R"| 7fc • 2; < 
?/, , V/c G /Cx} and P = {x £ R"! "yk ■ x <rij., V/c G JCx}- Then, the condition P^Q P (Z P translates 
to 7] < rjk < rif^, for k G Kx- Secondly, we will search the controller /i in a subspace spanned by a 
polynomial matrix: 

h{x) = H{x)e 
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where ^ S M"? is a parameter to be determined and the matrix H : M" — )> M^p^"?) is a given multivariate 
polynomial map. The use of a template expression is natural when searching for a controller with a 
particular structure. The input constraint (i.e. for all x S Rx, h{x) € U) is then equivalent to 

(2.3) VA: G /Cc/, Vx G Rx, au,k ■ H{x)9 < pu,k- 



Under these assumptions, the dynamics of the controlled system (12. 2p can be rewritten under the 
form 

±{t) = f{x{t),d{t)) + G{x{t),d{t))9, d{t) e D, 

where the matrix of polynomials G{x,d) = g{x,d)H{x). From the standard characterization of 
invariant sets (see |Aub91j ). it follows that the controlled system (]2.2p is P- invariant if and only if 

(2.4) VA; G JCx, Vx G Ffe,V(i G -D, 7fc • {f{x,d) + G{x,d)e) < 0. 

Then, Problem 12.21 can be solved by computing vectors 6 gM.'^ and rj G M'^' with r] < rjf^ < % for 
all k G fCx, and such that (j2.3p and (j2.4p hold. In the following, we first show how, given the vector 
rj G M' -^1 (and hence the polytope P), we can compute, using linear programming, the parameter 
6 (and hence the controller h) such that the controlled system ()2.2p is P-invariant. Then, we show 
how to compute jointly the controller h and the polytope P using an iterative approach based on 
sensitivity analysis of linear programs. Before that, we shall review some recent results on linear 
relaxations for optimization of polynomials over bounded polytopes |BG10] . 



3. Optimization of Polynomials over Polytopes 

In this section, we review some recent results of [BGIO] that will be useful for solving Problem 12.21 
Let us consider the following optimization problem involving a polynomial on a bounded polytope: 



minimize 


c • p{y) 




over 


y£R, 




subject to 


ai-y < bi, 


iel, 




aj -y = bj, 


J <^J, 



(3.1) 



where p : M™ — ;> M" is a multivariate polynomial map, c G M", R = [yi,yT] x • • • x [yrmUm] is a 
rectangle of M™; / and J are finite sets of indices; ak G M"*, bk G M, for all /c G I U J. Let us 
remark that even though the polytope defined by the constraints indexed by / and J is unbounded 
in R"^, the fact that we consider y £ R which is a bounded rectangle of M™ results in an optimization 
problem on a bounded (not necessarily full dimensional) polytope of M"'. Let p* denote the optimal 
value of problem (j3.ip . The approach presented in [BGlOj allows us to compute a guaranteed lower 
bound d* of p* . The approach is as follows. First, using the so-called blossoming principle |Ram89| . 
we transform problem ()3.ip into an equivalent optimization problem involving a multi-affine function 
on a polytope. The dual of this problem is then a linear program easily solvable and whose optimal 
value is a guaranteed lower bound of p*. 

3.1. Blossoming principle. Multi-affine functions form a particular class of multivariate polyno- 
mials. Essentially, a multi-affine function is a function which is affine in each of its variables when 
the other variables are regarded as constant. 
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Definition 3.1. A multi-affine function q : W^^ — )• M is a multivariate polynomial in the variables 
zi, . . . , zm where the degree of g in each of its variables is at most 1 : 



/(z) = q{zi, ..., zm) = Y^ qidu...,dM)4' 



, z 



Am 
M 



[M A „,„^ „ . TU,M 



where q(di,...,dni) ^ I^ for all {di, . . . , (ij\/) € {0, 1} . A map q : M — )■ M" is a multi-affine map if 
each of its components is a multi-affine function. 

It is shown in |BH06j that a multi-affine function q is uniquely determined by its values at the vertices 
of a rectangle R' of M^^. Moreover, for all x € R', q{x) is a convex combination of the values at the 
vertices so that we have the following result: 

Lemma 3.2. Let q : M*^ — )• M &e a multi-affine function and R' a rectangle ofM.^^ with set of vertices 
V' , then min (/(x) = vnvaqiv). 

The blossoming principle (see e.g. |Ram89j ) consists in mapping the set of polynomial maps to the 
set of multi-affine maps as follows. Let p : M™ — )• ffi" be a polynomial map. Let 6i, ... ,6m denote 
the degree of p in the variables yi, . . . ,ym respectively. Let A = {0, . . . , Ji} x • • • x {0, . . . , 5m}, then 
for all y S M"*, p{y) = {pi{y), ■ ■ ■ ,Pn{y)), where for all j = 1, . . . ,n, the components pj : M™' — )■ M 
are multivariable polynomial functions that can be written under the form: 

Pj{y)=Pj{yi,---,ym)= Yl Pj,{dr,...4rr.)yt ■■■v't 

{di,...,dm)eA 



+---+<5„ 



c<5iH 1-<5„ 



where Pj^(^d^,...,d„,) ^ I^> for ah (di, . . . , dm) G A and j = 1, . . . , n. 

Definition 3.3. The blossom of the polynomial map p : M™ — > M"" is the map q : M^ 
whose components are given for j = 1, . . . , n and z = {zn, . . . , zi Sj^, ■ ■ ■ ,Zmi, ■ ■ ■ ,Zm5 ) G 

by 

9i(^) = Yl Pj,(dl,...,dm)Bdi,5A^hl^ • • • ' ^1A) • • • Bd^,5miZm,l, ■■■, Zm,5rr.) 

{di,...,d,„)eA 
with 

Bd,5{zi,- ■-,25) = jjT Y/ ^<^i ■■■ ^'^d 

^d.) crGC(d,<5) 

where C{d, 5) denotes the set of combinations of d elements in {1, . . . , 5}. 

Example 3.4. The blossom of the polynomial map p{xi,X2) = {xi + xf, a^fx^) is q^ : M^ — )• M^ whose 
components are given by 

91(^1,1,2^1,2,^2,1,2:2,2,^2,3) = U^l,l + Zl,2) + Z2,lZ2,2Z2,3 

92(21,1, ^1,2, 22,1,^2,2, 22,3) = 2i,iZi,2|(22,lZ2,2 + 22,222,3 + 22,322,l) 

From Definition 13.31 it follows that the blossom q of the polynomial map p satisfies the following 
properties |Ram89j : 

(1) It is a multi-affine map. 

(2) It satisfies the diagonal property: q{zi, . . . , zi, . . . , Zm, • • • , Zm) = p{zi, . . . , Zm)- 
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(3) Let z,z' e M^i+-+^™, with z = (zi^i, . . . , z^^^, . . . , z^^i, . . . , z^^^^) and z' = (4^^, . . . , z^^^^, 
. . . , z^ ]^, . . . , z^ ^ ), we denote z = z' if, for all j = 1, ... , m, there exists a permutation 
TTj such that (2;j,i, . . . , 2j,<5j) = T^jiZjiT ■ ■ ^z', g,). It is easy to see that = is an equivalence 
relation. Moreover, for all z = z' , q{z) = q{z'). 

The diagonal property clearly allows us to recast problem ()3.ip for a multivariate polynomial map p 
as a problem involving its blossom q subject to inequality and equality constraints: 

minimize c • q{z) 
over z £ R', 

(3.2) subject to a/ • z < 6j, i € I 

a'j-z = bj, j € J, 

Zk,i = Zk,i+i, k = 1,... ,n, I = 1,... ,6j - 1. 

where R' = [yi,yi]^^ x ■ ■ ■ x [ym,ym]^"^ and the vectors a'^ are given for /c € / U J by 

' — ( ^'^'^ Qfc.i Qfc.m Qfc.m x 

3.2. Linear programming relaxation. In jBGlOj . based on Lemma 13. 2^ it is shown that the 
Lagrangian dual of problem (13. 2p is actually a linear program. After some complexity reduction 
enabled by the properties of blossoms stated above, the following result can be stated: 

Theorem 3.5 ( |BG10j ). The dual of problem l[3. 2\) is equivalent to the following linear program: 

maximize t 

over t £R, X € RI^I, /U € MI"^I 

(3.3) subject to Ai > 0, i e I 

t <c- q(y) + Y^ Xi{a[ ■v-bi)+Y^ fij{a'j ■ v - bj), v eV'. 
iei jeJ 

where V = V' / = with V' = {yi^yi} x • • • x {ym,ym} ■ Moreover the optimal value d* of this 
linear program is a guaranteed lower bound of the optimal value p* of problem 113. 1\} . 

The previous theorem provides a simple and efficient way to compute a guaranteed lower bound of 
the minimal value of a polynomial on a bounded polytope. In the following section, we will show 
how this can be used to solve our controller synthesis problem. 

4. Controller Synthesis 

In this section, given polytope P = {x € M"| Jk ■ x < r]k, VA; € JCx}, we show how to design a con- 
troller h : Rx — >• U such that the system (j2.2p is P-invariant. As explained before, we search 
the controller in a subspace spanned by a polynomial matrix: h{x) = H{x)9 where G M'^. 
Let F{x,d) = f{x,d) + G{x,d)9, it is a polynomial of degree 6i, . . . ,5n, pi, ■ ■ ■ , Pm in the vari- 
ables xi, . . . ,Xn,di, . . . ,dm- Its blossom is -Ffe = /fe + Gb9 where fb and Gb are the blossoms 
of / and G regarded as polynomials of degree 6i, . . . ,6n, pi, ■ ■ ■ , Pm- Let Hb be the blossom of 
the matrix H regarded as a polynomial map of degree 5i, . . . ,6n in the variables xi, . . . ,Xn- Let 

R^ = [Xi, W]^^ X-^[Xn, x:;^]^" , V^ = {Xi, XlY^ X • • • X {Xn, X^T}"^" , Vfj = {di, 7h}P^ X • • • X {drn, d^}^™ , 

Vx = Vy = and V[, = V'^l -. 

We first establish sufficient conditions such that h{x) E \J for all x € Rx'- 
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Lemma 4.1. If for all I € JCu, for allvx € Vi^, au,rHbivx)0 < Pu,h then for all x G Rx, h{x) S U . 

Proof. If for all I G K-u^ for all Ux G V^, at/,/ ' Hbivx)0 < /3;7,z, then using the third property of the 
blossom we have for all / G ICu, for all vx G V^, ajj^i ■ Hh{vx)d < f3u,i- Since ff;, is a multi-afhne 
map, Lemma [32] implies that for all I G /C^^, for all z G .R^, au,i ■ Hi,{z)9 < I3u^i. Then, using the 
diagonal property of the blossom, we obtain for all I G )Cu, for all x G Rx, ceu,i ' H{x)9 < I3u,i which 
is equivalent to say that for all x G Rx, h{x) € U. D 

The previous result gives a finite set of linear constraints which must be satisfied by parameter 9. 
We now establish conditions ensuring that the polytope P is an invariant for system (j2.2p . Let 
k G ICx, we will say that facet F^ of the polytope P is blocked if for all x £ F/., for all d £ D, 
Ik ■ if{^,d) + G{x,d)9) < 0. It is clear that system (|2.2p is P-invariant if and only if all facets are 
blocked. 

Lemma 4.2. Let 9 £ W^ and k G ICx, then the facet F^ is blocked if and only if the optimal value 
Pki^) ^f ^^^ following optimization problem is non negative: 

minimize —jk • ifi^, d) + G{x, d)9) 
over X G Rx,d G Rd 

(4.1) subject to aD,j • d < (3d., j G JCd, 

li-x < rji, i£lCx\ {k}, 

lk-x = rjk. 

A guaranteed lower bound d^{9) ofp'^{9) is given by the optimal value of the following linear program: 

maximize t 

over t gM, A'' gMI'^-^I,A'= gMI^oI^ 

subject to Xi > 0, i e ICx\ {k}, 

(4.2) ^ > 0. i e ^D, 
t < -7fe • ifbivx,VD) + Gh{vx,vT))9) 

\K.x\ \^d\ 

i=i j=i 



where for all i G JCx and all j G K-d vectors 7/ and a'^ • are given by: 

, _ /TM Till JbH 7iil\ ' — ('^R:M. "J^j.i ^D,hm Qp 

Ol Ol On On Pi Pi Pm P 



Proof. Remarking that Ff^ = Rx D F}^ and D = Ru n D, the first part of the Proposition is obvious. 
For the second part, let us remark that from the definition of the equivalence relation =, we have 
{Vx X V[))/ = that is the same as Vx' x Vd' ■ Then, we have just to apply the approach described in 
Section [3] where y = (x, d) and the multivariate polynomial p{y) is equal to /(x, d) + G{x, d)9. D 

Now we show how to choose 9 £ M3 such that the associated controller satisfy for all x G Rx, 
h{x) G U and the controlled system (12. 2 p is P-invariant. 
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Proposition 4.3. Let d* and (t*,{X'^*)k£K.x^i^''*)keKD^^*) ^^ ^^^ optimal value and an optimal 

solution of the following linear program: 

(4.3) 

maximize t 

over tGM, A^ gMI^^I,A^ gMI^oI,^gM^ k e ICx 

subject to Aj > 0, A; € JCx, i S JCx \ {k}, 

X';>0, kelCx, js/Cd, 

au,i ■ Hb(vx)0 < f3u,i, I e ICu, W G Vx^ 

t < -Ik ■ {fbivx,VD) + Gb(vx,W))0) 

\ICx\ \ICd\ 

i=l j=l 

Then, if d* is positive, the controller h{x) = H{x)9* satisfy for all x G Rx, h{x) G U and the 
controlled system i2.^) is P-invariant. 

Proof. We first start by remarking that problem (I4.3P is equivalent to the following optimization 
problem: 

maximize t 

over teR, 9 £R'^ 

(4 4) — 

subject to au,i ■ Hi,{vx)0 < I3u,i, I G JCu, vx G V^, 

t < dl{9), keJCx 

where d^(0) is the optimal value of linear program (|4.2p . Then, if d* > 0, this means that for the 
optimal 9*, we have for all k G ICx, dl*i^{9*) > 0. Therefore, by Lemma [4.21 all facets of P are blocked 
and thus the controlled system (j2.2p is P-invariant. The constraints on 9 also ensures, by Lemma [4.1[ 
that for all x € Rx, h{x) eU. D 

5. Joint Synthesis of the Controller and the Invariant 

In this section, we present an iterative approach for synthesizing jointly the controller h and the 
invariant polytope P solving Problem 12.21 It is based on sensitivity analysis of linear programs. At 
each iteration, we use a guess for the invariant polytope P. Following the approach described in the 
previous section, we try to synthesize a controller h that renders P invariant. If P cannot be made 
invariant by this approach, we use sensitivity analysis of linear program ()4.3p to modify P and obtain 
a new guess for the invariant polytope. The procedure is repeated until Problem 12.21 is solved. 

5.1. Sensitivity analysis. Let r/, // G M''^^!, let polytopes P = {x G M"] 7fc • x < r]k, Vfc G ICx} and 
P/j. = {a^ € ]R"| 7fc • X < r/fc + //fc, VA; G ICx}] P^ can be seen as a perturbation of polytope P. The 
main result on sensitivity analysis is given by the following proposition: 

Proposition 5.1. Let d* and it*, {X^*)keKx-' 0^^*)k(^K.D^ ^* ) denote the optimal value and an optimal 

solution of linear program ( [^.3[ ), let d* denote the optimal value of linear program |^.3| j where r] has 

been replaced by rj + ^, then 

d;>mmid*-X>^*.l,). 
keKx 
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Proof. For all k S ICx, for all vx G Vx' and V£) £ Vd', we have: 

-Ik ■ {fh{vx,vJ5) + Gh{vx,vJ5)0) + ^ Af*(7/ ■ vx - m - ^J^■^) + ^^ ^j*("Dj ' "^ - I^d,j) 

= -Ik ■ {fb{vx,vB) + Gb{v^,v5)0) + Y^ \t{li ■W-'ni) + Y^ ^^*(«D,i • vB - Pd,j) - \^* ■ ii 
> i*-A'=*-/i> mill (t*-A'''*-/x). 

k'&Kx 

This shows that ( min (t* — A * ■ ji), (A *)k^jCx^ (^ *)keK.D^^* ) i^ a feasible solution for linear pro- 
gram (|4.3p where 77 has been replaced by r/ + /U. It follows that d* > min (t* — X * ■ 12) which leads 

kSilCx 

to the expected inequality since d* = t*. D 

The previous result has the following implications. Let us assume that we are not able to synthesize 
a controller rendering polytope P invariant by solving the linear program (14. 3p . this means that 
d* < 0. Then, the previous result tells us how to obtain a modified polytope P^ in order to get 
dn>0 (or at least to get an improved d* > d*). This suggests that we can solve Problem 12.21 using 
an iterative approach described in the following paragraph. 

5.2. Iterative approach. Initially, let us assume that we have an initial guess for the polytope P; 
one can for instance use P but other choices are possible. We use an iterative approach to solve 
Problem I2.2| each iteration consists of two main steps. 

5.2.1. First step: synthesize a controller. Given polytope P = {x £ M^\ 'jk • x < i]k, \/k G ICx}, we 
use Proposition H3] to synthesize a controller h. Let d* and (t*,{X'^*)k£K:xj{^''*)keK.D^^*) denote the 
optimal value and an optimal solution of linear program (|4.3|) . If d* > 0, then we found a controller 
rendering P invariant for the controlled system ()2.2p and Problem 12.21 is solved. If d* < 0, then we 
move to the second step. 

5.2.2. Second step: modify the polytope. We now try to find /i € M' -^1 ensuring that polytope 
Pfi = {x £ M"'! 'Jk ■ X < r]k + fJ-k , yk £ ICx} will be invariant for the controlled system (j2.2p . For that 
purpose. Proposition 15.11 tells us that it is sufficient that d* — X^* • /-f > 0, for all k £ ICx. The 
requirement that P_Q P^ Q P can be translated to r]k — T]k ^ fJ'k ^ ^k — Vk for all k £ ICx- Also, 
since sensitivity analysis is pertinent mainly for small perturbations, we impose that for all k £ ICx, 
— £ < A*fc < £ where e is a parameter that can be tuned. Then, finding a suitable /i can be done by 
solving the following linear program: 

maximize t 
, . over t GM, /i G M'^-^I, 

^ ' ' subject to t < d* - X'^* ■ fi, k£}Cx, 

min(-e,?/^-r?fc) < Hk < max(-e, ?/fc - r?fc), k £ ICx- 

Let {t* , fi*) be a solution of this linear program. If the optimal value t* of this problem is non- 
negative then it is sufficient to prove that the controller h : Rx -^ U synthesized in the first step and 
polytope P^* solve Problem 12. 2[ Otherwise, if t* < 0, then we go back to the first step and start a 
new iteration with P = P^* . 
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Remark 5.2. Let us remark that the polytope P^. computed by solving (15. ip may have empty facets. 
In order to avoid such situations, it may be useful to replace fi* by fl* such that P^* has no empty 
facet and P^-* = Pji* (see Figure [T]). This can be done by solving a set of linear programs. 






Figure 1. The polytope P^* may have empty facets (center polytope), we replace 
//* by Jx* such that P^* has no empty facet and P^* = P^* (right polytope). 



Let us discuss briefly the computational complexity of our approach. Each iteration mainly consists 
in solving two linear programs. The linear program ()4.3p has 1 + q + \ICx\{\}Cx\ + \K^d\) variables 
and |/Cx|(|A^x| + |^Z)| + \^x\ + \^d\ ~ 1) + l^c/ll^l inequality constraints. Let us remark that 
|V^| = (5i + 1) X • • • X ((5„ + 1) and |y^| = (pi + 1) x • • • x (/9„ + 1). Since the complexity of linear 
programming is polynomial in average in the number of variables and constraints. It follows that 
the first step of the iteration has polynomial cost in the number of constraints of polytopes P, D 
and U and in the degrees of the polynomials. The linear program (j5.ip has 1 + |/Cx| variables and 
2|/Cx| inequality constraints. It follows that the second step of the iteration has polynomial cost in 
the number of constraints of polytope P. 



6. Examples 

Our approach was implemented in Matlab; in the following, we show its application to a set of 
examples. 



(6.1) 



6.1. Moore-Greitzer jet engine model. We tested our approach on the following polynomial 
system, corresponding to a Moore-Greitzer model of a jet engine |KKK95| : 

±2 = U. 

We first work in the rectangle Rx = [-0.2, 0.2]^ with disturbance d G D = Rn = [-0.02,0.02]. We 
want to synthesize a linear controller i.e h{xi,X2) = Oixi + 61X2 such that h{x) € U = [—0.35,0.35], 
for all X € Rx- Let P and P be polytopes with m = 24 facets with uniformly distributed orientations 
and tangent to the circles of center (0, 0) and of radius 0.01 and 0.2, respectively. Using our approach, 
we found the controller h{xi,X2) = 0.8076xi — 0.9424x2 rendering the polytope P shown on the left 
part of Figure [2] invariant. We make a second experiment, working in the rectangle Rx = [—0.2, 0.2]^ 
with disturbance d ^ D = Ru = [—0.025, 0.025]. We want to synthesize a polynomial control whose 
degrees are 3 in xi and 1 in X2 (i.e. the same as the vector field). P and P are polytopes with 
771 = 8 facets with uniformly distributed orientations and tangent to the circles of center (0, 0) and of 
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radius 0.01 and 0.2, respectively. Using our approach, we found a controller rendering the polytope 
P shown on the right part of Figure [2] invariant. The previous experiments show that by looking for 
controller of higher degrees, we may be able to find simpler invariants for larger disturbances. 



0.2 



-\/ /'///' -^ ^ '^ ■^ ■" 
0.15/ ////// ^ '■ ■^ 





Figure 2. Left: Invariant polytope P with 24 facets and a trajectory of (j6.ip illus- 
trating the invariance for disturbance d(t^ = 0.02cos(0.5t). Right: Invariant polytope 
P with 8 facets and a trajectory of ()6.ip illustrating the invariance for disturbance 
d{t) = 0.025 cos(O.U). 



6.2. Unicycle model. We now consider a simple model of a unicycle: 

X = vcos{ip), 
y = vs'm{ip), 

if = UJ. 

where v and w are the inputs of the system representing respectively the velocity and the angular 
velocity of the particle. In the following, we shall consider t; as a disturbance and co as the control 
input. Using the change of coordinates zi = xcos{(p) + ysin((/j), and Z2 = xsm(ip) — ycos{ip), we 
obtain the following polynomial system. 



■2) 



Z2 



V - Z2UJ, 
ZiiO. 



We work in the rectangle Rx = [-0.1,0.1] x [0.9,1.1] with disturbance v e D = Rd = [0.96,1.04]. 
We want to synthesize an affine controller h{zi,Z2) = Oq + Oizi + 02-22- In this example, we do not 
impose constraints on the value of the input. P and P are defined as polytopes with tti = 24 facets 
with uniformly distributed orientations and tangent to the circles of center (0, 1) and radius 0.01 and 
0.1 respectively. Using our approach, we found the controller h{zi^ Z2) = 1.0178 + 1.8721zi — 0.0253z2 
rendering the polytope shown on Figure [3] invariant. 

6.3. Rigid body motion. The last example is a model describing the motion of a rigid body. It is 
borrowed from |BI89] : 

{Xl = Ul, 
X2 = U2, 
±3 = X1X2. 

We work in the rectangle Rx = [-0.2,0.4] x [-0.2,0.2] x [-0.2,0.4]. In this example, we do not 
consider disturbances. We want to synthesize a multi-affine controller i.e /i : M^ — > M^ (defined by 
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Figure 3. Invariant polytope P with 24 facets and a trajectory of (j6.2p illustrating 
the invariance for disturbance v{t) = 1 + 0.04cos(0.1t). 

sixteen parameters), such that h{x) £ U = [—1, 1]^, for all x € Rx- Using our approach, we found a 
controller rendering the polytope with 18 facets, shown on Figure HJ invariant. 




Figure 4. Invariant polytope P with 18 facets and trajectories of (|6.3p illustrating 
the invariance. 



7. Conclusion 



In this paper, we have considered the problem of synthesizing controllers ensuring robust invariance 
of polynomial dynamical systems. Using the recent results of |BG10j on polynomial optimization 
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over bounded polytopes, we have developed an iterative approach to solve this problem. It is mainly 
based on linear programming and therefore it is effective. We have shown applications to several 
examples which shows the usefulness of the approach. Future work will focus on a deeper theoretical 
analysis of the properties of the linear programming relaxations of polynomial optimization problems 
as well as their application to other classes of problems in control. 
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